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ABSTRACT 



A goal of the NASA STEREO mission is to study the feasibility of forecast- 
ing the direction, arrival time and internal structure of solar coronal mass ejec- 
tions (CMEs) from a vantage point outside the Sun-Earth line. Through a case 
study, we discuss the arrival time calculation of interplanetary CMEs (ICMEs) 
in the ecliptic plane using data from STEREO/SECCHI at large elongations 
from the Sun in combination with different geometric assumptions about the 
ICME front shape (Fixed-^ (FP): a point and Harmonic Mean (HM): a circle). 
These forecasting techniques use single-spacecraft imaging data and are based 
on the assumptions of constant velocity and direction. We show that for the 
slow (350 km s" 1 ) ICME on 2009 February 13-18, observed at quadrature by 
the two STEREO spacecraft, the results for the arrival time given by the HM 
approximation are more accurate by 12 hours than those for FP in comparison 
to in situ observations of solar wind plasma and magnetic field parameters by 
STEREO/IMPACT/PLASTIC, and by 6 hours for the arrival time at Venus Ex- 
press (MAG). We propose that the improvement is directly related to the ICME 
front shape being more accurately described by HM for an ICME with a low 
inclination of its symmetry axis to the ecliptic. In this case the ICME has to 
be tracked to > 30° elongation to get arrival time errors < ± 5 hours. A newly 
derived formula for calculating arrival times with the HM method is also useful 
for a triangulation technique assuming the same geometry. 

Subject headings: Sun: coronal mass ejections (CMEs); Sun: heliosphere; Sun: 
solar-terrestrial relations; methods: analytical 
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Introduction 



Coronal mass ejections (CMEs) are episodic expul sions of huge amou nts of 



plasma and magnetic fields from the solar corona (e.g. 



Howard et al 



19851 ). CMEs 



Farrugia et al. 



1993 



Gosling 



1993 



Bothmer fc Zhukov 



e.g. 


Tsurutani et a . 


1988 


20061: 


Zhang et al. 


2007) 



Besides other detrimental effects, the resulting geomagnetic storms can knock out 
communication satellites, lead to errors in global positioning systems and induce cur- 
rents on the ground leading to wide-spread electricity blackouts (e.g. report by the 



Committee On The Societal &: Economic Impacts Of Severe Space Weather Events! 120081) . 



The NASA STEREO (Solar TErestrial RElations Observatory) mission ( 



Kaiser et al. 



20081 ) was launched in October 2006 and is designed to enhance our understanding of the 



initiation and propagation of CMEs, and thus foster our abil ity to forecast them , using 

20081 ). and 



simultaneous imaging in EUV and white-light coronagraphs (iHoward et al. 
sampling the solar w ind magnetic fi elds and plasma using the in situ instruments PLASTIC 



( iGalvin et al. 



20081 ) and IMPACT (ILuhmann et all 120081 ) . Additional to the novel three 



viewpoints through in situ instruments and coronagraphs (including SOHO / LASCO and 



the Wind/ ACE in situ measure ments near 



with Heliospheric imagers (His, 



Eyles et al. 



arth), both STEREO spacecraft are equipped 



20091 ) which can directly map the complete 



evolution of Earth-directed transients (such as ICMEs) in the inner heliosphere (< 1 AU) 
from a vantage point away from the Sun-Earth line. By definition, we use the term ICME 
for the structure after it has left any coronagraph's field of view (FoV), i.e. we use it for 
both the in situ and HI observations. 



Especially useful is the following concept based on earlier work by 



(119991 ) and demonstrated with STEREO/HI and SMEI by e.g. 



Sheelev et al. 



Rouillard et al. 



(2008 
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2009|), iDavis et al.l ( I2009h . Mostl et all (I2009bl l2010h and iHoward & Tappin! (I2009bh : if 
an ICME propagates with constant velocity and direction to large angles with respect to 
the Sun, an observer perceives deceptive angular accelerations/decelerations depending on 
its propagation direction relative to the observer. This effect can be used to assess the 
direction, velocity and arriv al time of an ICME by fitting known theoretical profiles to 



real-time ( IDavis et al.l 1201 ll ) or level 2 HI or SMEI data of the ICME's el ongation. There 



exist two practical realizations of this concept: (1) the original technique (jRouillard et al. 



20081 ) which assumes a point-like or negligibly narrow shape of the CME front (also callec 



the Fi xed-$ approximation; we call this Fixed-$ fitting or FPF), or the method by 



Lugaz 



( 120101 ) which assumes a circular CME front, attached to the Sun at all times at one 



end (Harmonic Mean Fitting or HMF). It has to be pointed out th at there also exist 



more co mplex approaches to the same problem like the T-H model (IHoward fc Tappin 



2009a 



b; 



Tappin fc Howard 



20091 ) where a large set of pre-existing solutions is fitted to the 



observations, also usefu l for real time application, and the white-light rendering method by 



Wood fc Howard! (120091 ) . In this study we use the FPF and HMF methods, and we denote 
the formulas for converting elongation to distance with FP/HM and the inverted method 
useful for predicting ICME parameters with FPF/HMF. 

The big advantage of the fitting techniques is that they are simple and fast to apply in 
real time, and they can be used for single spacecraft HI observations, which will most likely 
be available on future solar mi ssions or an observatory at the L5 point in the Sun -Earth 



system, 60° behind the Earth ( iSchmidt fc Bothmer 



1996 



Gopalswamy et al. 



201 lh . The 



big disadvantage is that these methods assume constant speed V and direction <p. In 
general, fast and more geoeffective ICMEs are expected to decelerate because the solar 
wind they run into is usually slower and exerts a drag force, and ICME di rections can be 



influenced by high speed solar wind streams (e.g. 



Gopalswamy et al 



20091 ). Thus, for a 



real-time application, the resulting speeds and directions will only be averages of the part 
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where the elongation profiles are measured, and their errors have to be carefully assessed. 
Using tw o HI observatorie s, the assumptions of constant velocity and direction can be 



dropped ( ILiu et al. 



2010aU b: 



Lugaz et al 



20101 ). and again there exist two limiti ng cases fo r 



the "e xtreme" geometric assumptions of a point or a cir cle: triangulation (TR 



Liu et al. 



2010al ) and tangent-to -a-sphere 



techniques is given by 



Liu et al. 



( TAS; 



feoiObh . 



mgaz et al 



2010l ). An excellent summary of these 



However, what is missing so far is a way how to calculate the speeds and arrival times 
for the methods assuming a circular front shape (HMF and TAS). It is the aim of this 
paper to derive an analytical formula which gives the arrival time and velocity of an ICME 
using the HMF and TAS methods. The best way to independently assess the validity of 
any of these techniques is to check the results against in situ magnetic field and plasma 
measurements of ICMEs. Thus, additionally, we assess the validity of our new formula by 
comparing the results of the FPF and HMF methods for an ICME e v ent on 2009 Febru ary 



13-18 (see also 



Patsourakos fc Vourlidas 



2009 



Kienreich et al. 



2009 



Cohen et al 



2009 



on the EUV wave), which we believe to be extremely well suited for this purpose for two 
reasons: the CME had a (1) practically constant velocity all the way in the HI FoVs and 
(2) it was propagating into a "simple" slow background solar wind, so that any differences 
in the results by the FPF / HMF methods should be based only on the assumption of the 
different geometry rather than other effects. 

T he pap er is structured as follows: in Section 2, we present an alternative version of the 



Lugazl (120101 ) formula for inverting the Harmonic Mean approximation for the global ICME 
leading edge shape, derived in Appendix A. Also, an associated correction formula for 
calculating the arrival time with HMF is discussed in Section 2 and derived in Appendix B. 
In Section 3, we discuss STEREO-A/El observations of the 2009 February 13-18 ICME 
event and application of the HMF and FPF methods. We check these results against the 
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ICME in situ observations by STEREO-B and Venus Express in Section 4. We summarize 
our conclusions in Section 5. 



2. Methods for fitting the time-elongation profiles of ICMEs 



2.1. Fixed-$ and Harmonic Mean 



For completeness, we revisit in this section methods to calculate a constant ICME 
direction and velocity from the observed elongation variation with time e(t) of a feature 
(o ften the ICME lead ing edge) in the H I field of views. For two excellent in-depth summaries 



sec 



Lugazl (120101 ) and 



LiuetaL 



(J2010b|) . We will restrict ourselves here to a short discussion 
and introduce definitions to be used in the following case study. We present in this section 
(1) a new alternative version of a formula which can be used for fitting time-elongation 
tracks of tra nsients assuming a circular front shape (Harmonic Mean Fitting or HMF, 
Lugazl (120101 )). and (2) a new geometrically consistent formula for deriving the ICME 
arrival time at a given location in the heliosphere with HMF. Both formulas are derived in 
the Appendix to this paper. 



The effects o 



by 



Sheeley et al. 



obse rving features of ICMEs at large angles to the Sun were introduced 



(11999T ). Basically, the assumption that the CMEs intensity signal originates 



from the plane-of-sky, as valid i n the coronaeraph fie. 



for larger angles from the Sun (jVourlidas &: Howard 



d of views close to the Sun, is n ot valid 



2006 



Howard & Tappin 



2009aj). The 

radial distance Rfp from the Sun of a point-like transient moving at a longitude <ft (<f) > 
is defined as westward) relative to the observer, who is situated at a distance d n from the 



Sun, which is observed at time t at an elongation angle e, is given by 
their equation Bl) as 



Sheeley et al. 



(1999 
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R FP {t) 



d Q sin e(t) 
sin(e(t) + 0) 



This is also known as the Fixed-0 (FP) approximation (IKahler fc Webbl 120071 ). These 
definitions are visualized in Figured! Solving for e(t) leads to 

Vppt sin( 4>fp) 



e{t) = arctan 
first written down explicitly by 



d — Vppt cos(4>fp) 



(2) 



Rouillard et al. 



( 120081 ). Here we have already replaced the 



distance Rfp by Vppt. This formula is thus only based on a constant velocity Vpp and 
a constant direction 4>fp- For comparison with real data, it is also necessary to obtain 
a launch time toFP on the Sun, for which e(toFp) = Oil- Equation <^ has the obvious 
advantage that a theoretical elongation profile can be fitted in a least-squares sense to 



the observed values 


o 


: e(t), which resu 


ts in the 3 


parameters 


/j? P ,d) W p, anc 


1 t()FP 


e.g. 


Rouillard et al. 


2008 


■ 


2009|; 


Davis et al. 


2009. 


2010; 


Mostl et al. 


2009b 




2010; 


Lueaz 


2010) 



We further call this the Fixed-(/>- Fitting or FPF approac h. We implemented the fitting by 



using a downhill simplex method (INelder fc Meadl (119651 ) . IDL function amoeba. pro) which 



is, from our experience, robust for this three-parameter minimization. It can also be done 



by creating an "error map" for Vfp and 



the theoretical and observed ICME track (j Williams et al. 



p and looking 



'or the minimum residual between 



2009 



Lugaz 



2010). 



The same pro cedure can be applied using the so-called Harmonic Mean approach. 



Lugaz et al 



(120091 ) derived a formula to convert elongation to distance, similar to 
Equation (j2J), but under the different assumption that the CME front can be approximated 
by a circle, which is attached to the Sun at all times at one end. This definition is 



l We note that for a consistent comparison of the launch time to solar surface signatures 
any e(t ) should be set to the angular size of the Sun's radius (~ 0.26° for an observer at 1 
AU), but to simplify the arrival time calculation e(io) — in this paper. 
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visualized in Figure 1. It is further assumed that the observer always measures the resulting 
elongation angle by looking alo ng the tangent to this circle (see Figure 1 and 



2009 



Howard & Tappin 



Lugaz et al. 



2009al ). The radial di stance of the circle 's apex Rhm, the point 



with largest distance from the Sun, is given by 



Lugaz et al. 



J2009J): 



RHM\t) 



2d a sine(t) 
1 ± xm(f(t) 4- (h)' 



(3) 



1 r 



An equivalent formula was also derived by iHoward fc Tappinl (j2009al their equation 41). 



We derived a new v e rsion of the inversion of this formula in Appendix A, compared to the 



one given by 



Lugazl (120101 ) . Setting Rhm = Vuhit and defining 
2d 



Vff ft/ft 



COS (pHM] 



sm{(p HM 



leads to 



,^ ,—b + ay/a 2 + b 2 - 1. 
e(t) = arccos( ^— ). 



(4) 



(5) 



This has only one term in contrast to the formula with two te rms gi v en by 



(120101 ). The reason we use this formula instead of Equation (1) by 



Lugaz! (120 lOf ) is that 



Lugaz 



we found it works more rob ust when run ning the downhill simplex method for the fitting, 



which we use in contrast to iLugazI (j2010l ) , who constructed error maps for a grid of Vhm 
and 4>hm values. It is our purpose to find out if the two approaches may produce different 
results in future studies and to further check the robustness of our fitting method, but in 



this paper we restrict ourselves to discussing the results of the simplex method alone. 



2.2. Arrival time calculation using fitting methods 

In this section we summarize the calculation of the arrival time t a of an ICME at a 
particular planet or spacecraft for both fitting methods. Beside the constant direction, the 
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FPF and HMF methods result in a constant velocity Vfp / Vhm and a launch time toFP / 
toHM on the solar surface, for which e(t ) = 0. For the FPF method the arrival time is then 
simply given by 



taFP — toFP + di/VFP, (6) 

with di being the radial distance of the in situ observing spacecraft from the Sun. For 
the FPF method, it is assumed that the ICME's white-light emission arises from a single 
plasma element, or in reality that the ICME front is of negligible small extent along the 
line-of-sight. But what does a "small extent" exactly mean? In this respect it is necessary 
to define for which cases an ICME with a direction 4>f p is actually expected t o hit a 



Davis et al 



2009 



201PJ) that if the 



spacecraft. It has been assumed in previous studies (e.g. 
direction cf>FP matches the in situ spacecraft longitude to within a few degrees, the ICME 
is expected to hit the spacecraft and the arrival time can be calculated by using Equation 
(J6]). This actually implies that the front shape of the ICME for the FPF method is a small 
segment (extending a few degrees in longitude) of a circle with radial distance of 1 AU from 
the Sun around the direction 4>fp- But for larger angular separations between the in situ 
spacecraft and the ICME direction this definition runs into troubles because it is unclear if 
the ICME is expected to hit the spacecraft at all. Clearly, only a statistical study using 
many ICMEs will determine where to set the boundary where a prediction should be issued 
that an ICME will hit a particular planet or spacecraft at a certain longitudinal separation 
angle (which we further call A) between the ICME direction and the target position. 

This problem is less evident for the HMF method because the ICME front shape is 
well defined as an expanding circle attached to the Sun at one end. However, the opposite 
problem arises here compared to the FPF method: there it is not clear if a spacecraft is 
hit at all by an ICME, but here sooner or later every in situ spacecraft should detect the 
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ICME front which is positioned at a longitude of 4>hm ± 90°, i.e., in the half space centered 
on <Phm- Nevertheless, the well-defined circular geometry directly leads to an analytical 
arrival time calculation, which we derive in Appendix B. To calculate an ICME arrival time 
at a given heliospheric position which is geometrically consistent with the circular front, it 
is necessary to correct for the fact that all parts of the circular front except for the apex 
will always arrive later at a given radial distance from the Sun than the apex (see also 
Figure [8^). To calculate the arrival time of an ICME at a given heliospheric distance di and 
at a longitudinal separation A with the HMF method, an arrival time correction t c must be 
added to the arrival time of the ICME apex: 

taHM = toHM + 77 H t C - (7) 

As derived in Appendix B, this can be reduced to the simple equation 

taHM = toHM + 77 7~TT • (8) 

V HM cos(A) 

It is seen in a straightforward way that the circular HM geometry leads to the elementary 
result that the speed of the ICME flank at a given angle A to the ICME apex is just 
reduced by a factor cos(A), for this particular geometry. This is especially important for 
comparison to in situ measurements of the solar wind bulk speed when using the HMF 
method, where the speed Vhmi = Vhm cos(A) has to be used, in contrast to FPF, where 
just Vpp has to be taken. 

Additionally, there exists a triangulation version of the HMF approac h, which was 
called the tangent-to-a-sphere (TAS) method (Model 1 in Lugaz et al. 2010 ). In contrast to 
the single-spacecraft fitting methods, using two e(t) functions gained by widely separated 
observer s makes it possib le to relax the assumptions of constant speed V and direction </> 



(see also 



Liu et al. 



2010a 



bj). For calculating the ICME arrival time with TAS, Equation ([8]) 



also has to be used, but because VpAsif) and 4>TAs{i) are now functions of time a final speed 
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and final direction have to be put into Equation (jSj), i.e. the resulting speed and direction 
values at the last time step both e(i) functions could be determined. We imagine that the 
final direction 4>tas is simply the result of the last e(t) pair, but for the terminal propagation 



speed Vtas averages of the last few hours need to be taken 



Decause of the fluctuations 



Liu et al. 



2010a|). Alternatively, the 



resulting from the derivatives (see e.g. Figure 4 in 
arrival time correction t c , which is given by Equation (1B5P in the Appendix B in easy-to-use 
units, needs to be added to an already calculated arrival time of the ICME apex at a given 
radial distance from the Sun. 

Before proceeding to a data example, we also note that we extract the ICME track 
from STEREO-A/B.I images in the same plane (the ecliptic) as the main in situ observing 
spacecraft (STEREO-B) is situated. In this way we avoid any effects based on observations 
and positions yielding from different latitudes. 



3. Case study of the 2009 February 13-18 ICME 



On 2009 February 13 a CME er upted from a small bip olar active region (AR 11012) at 



position S04E46 as seen from Earth (IKienreich et al. 



2009 



Miklenic et al. 2011, submitted 



to Solar Physics). It was associated with a GOES class B2 . 3 flare peaking at 05:46 UT and 



(2009 


) and 


. ~, ^ - 

Cohen et al. 


(2009) 



Kienreich et al 



( I2009h : 



Patsourakos &: Vourlidas 



average at 1.0033 AU and at a longitude of —47.6 degrees east of the Earth in Heliocentric 
Earth Ecliptic (HEE) coordinates, STEREO-A at 43.5 degrees west of the Earth at 
0.9643 AU, and their separation in HEE latitude was negligible (< 0.3°). This means that 
the event was seen at the limb from STEREO-A and almost exactly at disk center by 
STEREO-B, making this a practically perfect quadrature configuration, as the STEREO 
longitudinal separation was 91 degrees. The two probes are equipped with instruments 
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(SECCHI, iHoward et al.l 120081 ) to image the Sun's corona in the extreme ultraviolet 
spectrum (EUVI), and the inner (COR1, 1.5-4 R ) and outer corona (COR2, 2.5-15 R ) in 
white-light. Images of t he heliosphere, a 



Imagers (SECCHI/HI; 



Eyles et al. 



so in white-light, are obtained by the Heliospheric 
20091 ). with HI1 covering 4-24 degree elongation from 
the Sun and HI2 18.7-88.7 degrees. These values are valid for the ecliptic plane, at which 
the optical axes of the HI telescopes are pointing nominally. STEREO is also able to 
probe the solar wind at the pos ition of each spac ecraft itself with tw o in situ instruments: 



IMPACT flLuhmann et al 



20081 ) and PLASTIC flGalvin et al. 



2008|) 



Miklenic et al. (2011, submitted to Solar Physics) discussed the early acceleration 
phase of the event and found that the CME reached its propagation velocity of 350 km s _1 
already at around 2.5 R , which is still in the field-of-view (FoV) of the STEREO/COR1A 
coronagraph and thus rather close to the Sun and well below the HI FoVs. The ICME was 
detected in situ by STEREO-B and Venus Express (see Section 3.2), with an observed 
average proton bulk speed at STEREO-B also of ~ 350 km s _1 . The state of the 



background solar wind for this event is discussed by T emmer et al. 



'2011, submitted to 



20031 ). who found a rather 



ApJ) using the WSA/MAS+ENLIL solar wind model (lOdstrcill 
simple state of the heliosphere around the time and location of the eruption, with a small 
high-speed-stream to the east of the source region. STEREO-B EUVI images also show that 
the AR 11012 was the only one on the solar disk viewed from the location of STEREO-B. 

In summary, the CME was slow, reached its terminal velocity very close to the Sun and 
propagated into a simple, slow solar wind, approximately maintaining its constant velocity 
from 2.5 R up to 1 AU (215 R ) (see also Rollett et al., 2011, submitted to Solar Physics). 
These characteristics make the event especially well suitable for a comparison between 
the HMF and FPF methods, because the ICMEs propagation in the HI FoVs, where the 
methods are applied, will not be affected strongly by any interaction with fast streams from 
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coronal holes (e.g. 



Gopalswamy et al 



20091 ) or other ICMEs, thus minimizing the effects 
of a possible deflection or acceleration during its propagation through the STEREO-A 
FoVs, which would be in conflict with the FPF/HMF assumptions of constant velocity 
and constant direction. Thus any differences in the results may be safely ascribed to the 



strongly different assumptions on the shape of the ICME front as used for FPF/HMF. 



3.1. STEREO-A heliospheric imaging observations 



In Figure [2] the top 2 images show running differences of the ICMEs evolution in 
STEREO-H.11A and the bottom 2 images show HI2A. The intensities in the images 
are related to Thomson-scattering of white-light solar radiation off solar wind electrons, 
and the incident i ntensity for a given observer maximizes along the li ne-of-sight at the 



Thomson-surface (jVourlidas fc Howard 



2006 



Howard fc Tappin 



2009al . see also Figure |6]) . 



We encourage the reader to also look at the HI1A/HI2A movie accompanying the paper in 
the electronic edition of The Astrophysical Journal. Overall, the CME is rather small and 
faint. It is actually so faint that it cannot be discerned in STEREO-B HI movies, thus the 
observations limit our study to using single-spacecraft HI methods, so it is not possible to 
use any triangulation approaches. In the first image, it is seen that the ICME clearly has a 
ringlike shape with a small angular w idth of roughly 20°, which is consistent with a view 



of a croissant-shaped ICME edge-on (IThernisien et al 



2009 



Wood fc Howard 



2009J), i.e. 



we look along the axis of the ICME, which is pointing out of the paper plane. The second 
density enhancement at its sunward end disappears around —15° elongation, leaving only 
the leading edge visible. 



We estimated the ICMEs angular width (AH 7 ) in HI using the method of 



Lynch et al. 



( 120101 ) (results denoted by AWl) and additionally by a direct measurement in the 



plane-of-sky (POS) of STEREO-A (denoted by AWd). The results are indicated in Figure [2] 



-14- 



for each frame and show that between early Feb 14 and late Feb 15 AW stayed between 
19 and 27 degrees, which is practically constant. In the accompanying online movie it is 
also seen that neither the angular width nor the latitudinal propagation direction and the 
morphology significantly change during the full propagation in the HI FoVs. Later than 
early Feb 16 the ICME becomes very faint in a running difference image, but its leading 
edge still can be seen in the Jmap (see below and Figure Eb). 

In the right two images of Figure |2] it is also seen that the ICME propagates along a 
central position angle (PA) about 2 degrees north of the ecliptic plane, and the latter is 
drawn at a PA of 92.0° in the STEREO-A POS. The ICMEs central PA in the HI1 frame is 
87.5°, and 90.5° in the HI2 frame. Thus, in latitude, the ICME propagates outward close 
to the ecliptic, and thus close to STEREO-B (at PA = 92.0°) and Venus (at PA = 90.3°). 
The heliospheric longitude of STEREO-B is only 1° to the east of the Solar source region of 
the CME, and Venus is located 21° to the west of it, at a radial distance of 0.7185 AU and 
a HEE longitude of —25.4 degrees (roughly halfway between STEREO-B and the Earth). 
Thus, judging from the location of the Solar source region, the spacecraft positions and the 
symmetry axis close to the ecliptic plane the ICME is expected to hit both STEREO-B and 
Venus. 

Figure [3^ shows the proton density N p observed by STEREO-B. Panel (b) shows a 
Jmap extracted along the ecliptic plane. This is done by stacking vertically the median 
of the central row s ±16 pixels in the STEREO-A HI1/2 running difference images (for 



further details see 



Davies et al 



20091 ). It is seen that on early 2009 February 18 a rise in 
N p corresponds to the arrival of a track in the Jmap which is the leading edge of the ICME 
seen in Figure |2j The track in the Jmap was manually selected by choosing 30 points along 
the front of the leading edge, for 10 times. From the uncertainties associated with the 
manual selection of the points we will discuss how the results of the fitting procedures are 
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affected, so all systematic errors are assumed to arise from the manual selection only. The 
measured points were then interpolated onto an equidistant temporal grid consisting also 
consisting of 30 points (see Figure [3b) . 



3.1.1. Results of the fitting methods 



We now applied the FPF and HMF methods, as introduced in Section 2, to the 
measured ICME track. Figure [3b shows the results. In panel (c), the fitting functions (red: 
FPF; blue: HMF) are shown along with the observed e(t) and the error bars arising from 
the manual selection of the track. Both methods reproduce e(t) very well, but the results 
for the free parameters differ considerably (see also Table 1): 



(j) FP = _79°, v FP = 280 km s -1 , 
W = -107°, V HM = 322 km s" 1 , 



(9) 
(10) 



so the HMF direction is 28 degr ees eas t ward (or further away from the observer) compared 



to FPF. This has been noted by 



Lugaa (120101 ) to be a general feature - especially for a limb 



event the observer will mostly see the flank of the ICME in HI because of the scattering at 
the Thomson surface. By "flank" we denote those parts of the ICME front which are closer 
to the Sun compared to the apex, the point along the circle with the greatest distance from 
the Sun. The HMF method corrects this effect and will put the ICME apex direction always 



farther away from the observer compared to FPF (ILugaz 



2010j). In this case the directions 



for FPF is 12 degrees west of the longitude of STEREO-B and HMF 16 degrees east. 

The expected in situ speed at the location of STEREO-B of Vhmi — 310 km s _1 is 
30 km s _1 faster compared to FPF. We attribute this again to the fact that for a limb ICME 
event, especially one with low axis inclination to the ecliptic (to be justified further in the 
next section), more and more of the ICME flank is sampled. For a circular front, points 
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at the flank will always move slower than the ICME apex, and again the HMF techniqu e 



2011) 



corrects this effect which results in a higher speed compared to FPF ( ILugaz et al. 
Because the launch times to are quite similar (see Table 1 for a summary), the different 
speeds lead to quite different arrival times: due to its higher speed, t a HM, which already 
includes the arrival time correction we presented in Section 2, is about 12 hours earlier 
compared to FPF, which is a considerable difference. 

In Figure [3H we demonstrate how well both fits are able to reproduce the observe d 



function e(t). To this end, we use the fitting residual ([Williams et al. 



2009 



Lugaz 



2010) 



N 



a 2 = l/Nj2«U)-e flt (U)) 2 (11) 
i=i 

with N being the number of data points along the ICME track and eju the calculated 
elongation value for the time step U using the given set of fitting results V, and to- For 
the least-squares minimization using the downhill-simplex-method, a slightly modified a 2 
is used, where we replaced both e terms with their norms |e(tj)| and |e/it(*i)|> which lead 
to a more robust minimization than Equation ( fTTj) . From Figure [3H, where we plotted 
e(tj) — €fit(ti) for every e(t) measurement, it is seen that both functions are able to almost 
equally well describe the observed elongation variation of the ICME leading edge. The 
residual a 2 is a factor of 2 higher for FPF than for HMF, and both residual functions 
follow each other very closely, with the strongest deviations to the observed elongations 
at the end of the track. Nevertheless, based on the fits one cannot clearly distinguish 
between the two methods based on the consideration that a goodness-of-fit is better for 
one method as compared to the other. A real-time space weather predictor would thus not 
be able to distinguish if one of the predictions should be preferred over the other because 
he/she is able to fit the e(t) function with almost the same accuracy with both methods. 
To the contrary, both functions have three free parameters and are flexible enough to fit 
e(t) almost perfectly, but the results of the inversion differ quite greatly due to the vastly 
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different assumed geometries. Thus, the only way to distinguish between the success of the 
two methods is to look if the independent in situ observations, available for this event even 
by two spacecraft, are consistent with the results derived from the heliospheric images. But 
before we proceed to discuss the in situ data, we assess the errors involved in the derived 
parameters. 



3.1.2. Error analysis 

How do the errors arising from the manual selection of the ICME track in the Jmap 
impact the results of the FPF and HMF methods? To this end, we constructed two 
additional tracks 

e+(t)=e(f)+2<7 e (f) (12) 
e_(*) = e(t)-2<7 e (*) (13) 

with (T e {t) being the standard deviation of e(t). Thus e + (t) / e_(i) is the track with 
the maximum positive / negative deviation from e(t). We again fitted the two tracks of 
maximum deviation with the FPF and HMF methods and averaged the difference between 
these results and the ones for e(t) (Table 1) to obtain an error for V, <p, to and t a . For the 
latter the HMF arrival time correction has to be included. We found a mean standard 
deviation of 2a e = 0.51° arising from the manual selection of points, thus the systematic 
error associated with the manual measurement in the Jmap is about ±0.5°. This error 
affected the direction as <Pfp = —82° ± 1.8 and 4>hm = —111° ± 2.8. The errors in the 
speeds were Vfp = 288 ± 5.0 km/s and Vhm = 332 ± 8.7 km s, so again slightly higher for 
HMF. The error in launch times were ± lh 35min (FPF) and ± lh27 min (HMF). Most 
importantly, the arrival time errors are ± 54min (FPF) and ± 20 min (HMF), where we 
quoted the values for the location of STEREO-B. Interestingly, the error is less for the 
HMF method by about a factor 2, which arises from the property that the arrival time 
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correction partly compensates for the error in the arrival time due to the slightly different 
(f). It is important to note that the arrival time error of less than 1 hour is much shorter 
than the difference of 12 hours in t a between the two methods (see Table 1). It is clear that 
these errors are rather small because we used a very long track which reaches completely 
up to the elongation of the spacecraft in the Jmap. We will test the accuracies and errors 
of shorter tracks in Section 13.31 but before we check the results obtained in this section 
against independent multi-point in situ data. 



3.2. STEREO-B, VEX and ACE in situ observations 



The STEREO-B spacecraft was situated at -91 degrees (HEE longitude) to STEREO- A 
during our event. This is A# = —9° and A# = +20° away from the ICME directions given 
by FPF and HMF, respectively. Thus, the directions from both methods in dicate that 
CME will hit STEREO- B. Figure H] shows the proton bulk parameters (IGalvin et al 



the 



20081 ) and the magnetic field (jLuhmann et all 120081 1 measured in situ. Because of its low 



velocity, the ICME does not drive a clear shock, but between 2009 Feb 18 08:00-10:00 UT 
two shock-like features mark the transition from very slow solar wind (300 km s _1 ) to an 
enhanced density region traveling with 350 — 360 km s _1 . We thus define the arrival time 
at STEREO-B £g= Feb 18 10:00 UT, which marks the beginning of the proton density 
enhancement (first vertical guideline). Because we tracked the leading edge of the ICME 
in the Jmap, also related to the most enhanced densities, this is the time which has to 
be used for assessing the arrival time predictions. After what could be a small-scale flux 
rope a discontinuity around Feb 18 17:00 UT marks th e beginning of a larg e-scale flux rope 



(diameter 0.186 AU), which is not a magnetic cloud (IBurlaga et al.lll98ll ) because the 
magnetic field is not smooth and the temperature is higher compared to its environment. 
Nevertheless, visually it can be perceived that the field rotates as the component (in this 
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case the poloidal one) goes from north to south and Bt (the axial component) is negative 
or eastward, throughout the flux rope interval. Assuming this rotation to be caused by 
the spacecraft going through a helical flux rope, it is of north-south-east (NES) type and 
of right-handed chirality accor ding to the classification scheme for bipolar magnetic clouds 



by iBothmer fc Schwennl (1l998l ). We conclude that the in situ magnetic field and plasma 
parameters at STEREO-B are in rough agreement with what is seen in the HI images: a 
slow ICME with an axis inclination 9 close to the ecliptic, and hit slightly off-center because 
of its northward propagation direction. 

In Figure H] we also plotted the arrival times of the HMF method (vertical dashed blue 
line) with the velocity Vhmi (horizontal blue line), and the same has been done for the FPF 
results in red. The difference between the predicted arrival time by HMF and the actual 
arrival time at STEREO-B is t a HM ~ — +8 hours, so positive values indicate arrival 
times later than predicted. For, FPF, t a pp — t B = +20 hours, so the predicted arrival is 
almost a day later than the observed one. Here, the new HMF method including the arrival 
time correction would have given a much better prediction. The difference to the actual 
arrival time is also well outside of the error bars which arise from the manual selection 
of the ICME front. We thus attribute the residual discrepancy mainly due to the ICME 
leading edge geometry, which could clearly be different from the Harmonic Mean circle, 
and associated line-of-sight effects. Nevertheless, the point is that the new HMF method is 
more accurate by about 12 hours for this particular ICME than the original FPF method. 

Venus Express (VEX) was positioned —69° away from STEREO-A in longitude, so 
A v = +13° and A v = +42° away from the ICME directions given by FPF and HMF, 
respectively. While one will expect the ICME to hit VEX from the FPF method, with a 
direction between STEREO-B and VEX, it is not immediately clear if the HMF direction 
indicates a hit, though, as already discussed before, the HMF geometry, if taken literally, 
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will predict a hit at any spacecraft in the half space centered on 4>hm- The arrival time 
correction for Venus Express is already about 30 hours at this longitudinal separation. I n 



Figur e we show B measured by VEX in Venus solar orbital coordinates (jZhang et al. 
20061 ). but rotated by 180 degrees (B x = —B x yso]B y = —B y vso) so it approximately 



matches the RTN coordinates orientation we used for STEREO-B. We linearly interpolated 
over time intervals when VEX was below 3 Venus radii to see only the solar wind magnetic 
field. A magnetic flux rope interval, with a smooth rotation during an elevated field 
strength of about 10 nT, is found to be starting at 2009 Feb 17 22:00 UT, which we define 
as the Venus arrival time ty. The interval continues until about 2009 Feb 18 14:00 UT. 
Qualitatively, its field is oriented mostly southward and westward, but not clearly rotating 
and also not consistent with the field rotation at STEREO-B. However, we think it is very 
unlikely that this is connected to a different ICME event, because the ICME is seen passing 
over the location of Venus, which is very close to the location of STEREO-B along the 
line-of-sight, in the online movie around Feb 18. Also, in latitude, the ICME is expected to 
hit Venus centrally. This puzzling inconsistency in the field could point out that the part of 
the flux rope where the field is most strongly wound and which carries the major part of the 
magnetic flux and helicity away from the Sun does not extend over its full length back to 
th e Sun, but is rather confined around the apex part of the ICM E (recently also discussed 



by 



Yamamoto et al. 



2010 



Mostl et al 



2010 



Kahler et al. 



201lh . The error in the arrival 



times for FPF was t a pp — ty ~ —10 hours and for HMF t a uu ~ ty ~ +4 hours, which 
means that again the consistency with the HMF method is better. However, we note that 
without high time resolution data for solar wind plasma parameters it is rather difficult to 
set the arrival time unambigously and the arrival of the ICME is less well defined as for 
STEREO-B. 

In Figure [6] we plotted the event geometry and visualized the ICME directions and 
spacecraft positions. It can be seen that the circular front assumed by HMF leads to a 
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negative shift in £ a compared to FPF at the position of STEREO-B and to a positive shift 
in £ a at the position of VEX, as a direct consequence of the different geometries. Both time 
shifts from the FPF to the HMF method are more consistent with the in situ arrival times 
as discussed above. The improvement in arrival time consistency is 10 hours at the position 
of STEREO-B and 4 hours at VEX. It is seen in Figure [6] that this arises from the circular 
front which touches the STEREO-B and VEX positions, despite the very different radial 
distances, at almost the same time (see also Figure [5]). 

The Earth was —38.5° and —67° away from the ICME directions given by FPF and 
HMF, respectively, at a heliospheric distance of 0.9877 AU. While FPF gives an arrival 
time of 2009 Feb 18 22:53 UT at the ACE /Wind satellites, a value which is quite close to 
the arrival time at STEREO-B because of the almost similar heliospheric distance, the 
arrival time given by HMF is 7.5 days later on 2009 Feb 26 07:49 UT, which is due to a 
large arrival time correction because of the large longitudinal separation between Earth 
and the HM apex. We omit the plot showing the in situ data by the ACE /Wind spacecraft 
near Earth because there are no clear signatures of an ICME or any large-scale magnetic 
flux rope between Feb 18-26. This means that a longitudinal separation of 40-70 degrees 
is too far away to detect the flux rope in this case, even if it has a low axis inclination. If 
we assume the ICME apex at 1 AU is close to STEREO-B (so the ICME travels radially 
outward), and the ICME has an axis inclination to the ecliptic 9 « 0, then the non-detection 
near Earth 47° west of the apex leads to the constraint that the longitudinal width of 
this slow solar-minimum IC ME must be < 94°, which is roughly consistent with the value 



of « 70° (jWood et al. 



2010) for a different slow solar-minimum ICME event observed by 



STEREO /HI, and the < 60 ° extent obtained from multi-spacecraft in situ observations by 



Bothmer &; Schwennl ( 119981 ) 
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3.3. Accuracy as a function of track length 

So far we have focused on the systematic differences which arise from the different 
geometric assumptions underlying the two techniques. But a major goal of future missions 
concerned with real-time space weather forecasting should be to have the longest possible 



prediction lead times. Comm only it is quoted that the ICME must 



30 degrees from the Sun (e.g. 



Williams et al. 



2009 



De observed to at least 



Lugaz et al. 



20101 ) to get reliable results 



with HMF and FPF. To check this statement for our case study, we plotted in Figure [T^l 
the error in the arrival time at STEREO-B as a function of the angle of the maximum e 
measurement taken as an input for the fitting methods. We also include the error bars 
for t a described in Section 13.1.21 If only the measurements from the HI1 instrument are 
used between e = 4 — 24°, the resulting errors in the arrival time are ±12.5 hours, so 
basically a day, for FPF, but much less for HMF, due to the compensation by the arrival 
time correction mentioned above. As expected, for higher values of max(e) the errors, in 
particular for FPF become less and less, and are ±4.8 hours for max(e) = 30°, corresponding 
to a prediction lead time of about 48 hours until the actual arrival for this slow ICME. 
Also, for max(e) > 30°, the predictions by FPF are always delayed with respect to the HMF 
values by ~ 10 hours, thus, the systematic difference between the arrival times given by the 
two methods is not affected by the choice of max(e). 

In Figure [7b, we plotted in a similar manner the difference between the heliospheric 
longitude of STEREO-B and the propagation directions as results of the fitting methods for 
increasing values of max(e). We do not know the exact propagation direction of the ICME, 
but looking at the error bars, their size becomes clearly more reasonable for max(e) > 30°, 
where all directions range between A# = [—20, 10]°, and they would all have predicted a 
hit at the STEREO-B spacecraft. Also, for all times, the HMF method gives a direction 
further away from the observer than FPF, and the size of the error bars evolves similar in 
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time. 



4. Discussion and conclusions 

We discussed several new elements concerning the relationship between heliospheric 
images of interplanetary coronal mass ejections obtained by STEREO/HI and their 
signatures observed in situ. We subjected two methods to forecast the ICME parameters 
direction, speed and arrival time to close scrutiny. These methods approximate the ICME 
front by a point (Fixed-$) or a circle (Harmonic Mean), and the observed elongation 
of the ICME front as a function of time is fitted with inverted formulas to derive these 
parameters. We introduced (1) an alternative, simpler version for the inverted formula for 
the Harmonic Mean fitting method, and (2) an analytic arrival time correction formula 
which is geometrically consistent with an assumed circular front shape of the ICME. It 
should be applied to any real time predictions or case and statistical studies which discuss 
ICME arrival time calculation using met hods which assume a circular ICME front, i.e. 



single-spacecraft H armonic Mean Fitting (ILugaz 



(ILugaz et al. 



2010). 



20101 1 or two-spacecraft tangent-to-a-sphere 



We then checked the validity of the derived formulas by applying them in a case study 
to the 2009 February 13-18 ICME event, a perfect quadrature event, seen at the limb by 
STEREO-A and in situ by STEREO-B and Venus Express. The event was a slow event, 
with approximately constant speed throughout the HI field of views. We found that the 
arrival time given by the Harmonic Mean Fitting (HMF) method was more accurate than 
the Fixed-<3> fitting (FPF) method by 12 hours, in direct comparison to the observed in 
situ arrival time at STEREO-B. We attribute this to the fact that for an ICME with 
low axis inclination 6 to the ecliptic, the ICME front is better described by the circular 
geometry assumed by Harmonic Mean. The low 9 value was inferred independently by the 
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morphology and small de-projected angular width of the ICME in STEREO-A/HI and 
the magnetic field rotation observed by STEREO-B. At the latter spacecraft, the rough 
behavior of the magnetic field components pointed to a slight off-center hit of the ICME, 
consistent with a northward propagation by a few degrees of the ICME seen in HI. Thus, 
STEREO-B crossed the outer part of the cross-section of the flux rope. These observations 



are also consistent with a view that the outer part of the flu x rope is c 



and not a magnetic cloud, with a rather flat pressure profile ( IJian et al. 



assified as an ejecta 



2006|). Additionally, 



Venus Express observed a magnetic flux rope too, about 23 degrees west of STEREO-B, 
and again the in situ arrival time was more consistent with the one predicted by HMF. We 
also assessed the errors in the arrival time predictions as a function of how long the ICME 
is visible in the heliospheric images, and found that one needs to track the ICME up to 30 
degrees in elongation to obtain systematic errors (arising from the manual selection of the 
track only) which are less than ±5 hours when using scientific (level 2) HI data. 

It is our new hypothesis in this paper that the two methods for fitting time-elongation 
data (HMF and FPF) are just limiting cases of ICMEs with low or high axis inclination 



to the ecliptic, using two "extreme" assumptions on the ICME front sh ape (circu 



point-like). The same can also be said for the trian gulation (point-like 



Liu et al. 



ar vs. 



and tangent-to-a-sphere (circular, 



Lugaz et al 



2010af ) 



20101 ) methods, which use two-spacecraft 



HI observations. In our hypothesis, the HMF method corresponds to the case of low 
axis inclination and thus a large longi tudinal extent in the ecliptic (e.g. 70 degrees for a 



solar minimum low-inclination ICME, 



Wood et al. 



201(1 1 along the line-of-sight. The FPF 



method should give better results for an ICME with a high axis inclination because the 
extent of the cross-section is generally thought to be much smaller along the line of sight 
than its lo ngitudinal extension (e.g. about 20 degr ees for a solar minimum high-inclination 



ICME, see 



Kilpua et al 



2009 



Mostl et al. 



2009al lch. 
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This event was especially well suited for testing this hypothesis because the constant 
speed and constant direction assumptions of the HMF and FPF methods were approximately 
fulfilled, thus the systematic differences we found should arise only from the truly different 
shape of the ICME front. In reality, these simplified geometric assumptions will only be 
rough approximations the ICME front shape, and for slow events, which tend to get swept 
up in corotating interaction regions and become more easily distorted, any simple geometric 
model might be too restrictive in general. But for faster and more geo-effective events, 
where the front shape might be expected to be less easily affected by the background solar 
wind structure, the fitting methods could provide increasingly accurate results. 

We plan to extend this type of study for a larger set of events to find out how feasible 
these techniques are, if they have to be extended to other geometries of the ICME front and 
to further test our hypothesis concerning the differences in axis inclination. Also, one of the 
major questions wil l be if one or rather two spacecraft are nec essary for a L5 or a combined 



L4/L5 mission (e.g. 



Liu et al 



2010b 



Gopalswamy et al 



20111 ) to provide us with sufficient 



accuracy for advance warnings of solar eruptions using heliospheric imagers. The methods 

so for any upcom ing space 



Davis et al. 



(1201 lh showed 



described in this paper can be used for such a mission, and al 
weather forecasts at Earth and other planets with STEREO. 
the efficacy of using the FPF method for real-time prediction of an ICME event in April 
2010, when the STEREO spacecraft were separated by about 70 degrees to Earth, which is 
comparable to the L5 point at 60 degrees. There is still a long way to go, but with a solar 
maximum with many ICMEs in front of us we will be able to use STEREO extensively to 
continue honing our skills. 
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Table 1: A compilation of observed parameters and results of fitting techniques for the 
2009 February 13-18 ICME. The results for COR1A are taken from Miklenic et al. (2011, 
submitted to Solar Physics). FPF: Fixed-0 fitting; HMF: harmonic mean fitting; The suffix 
"/B" indicates results for STEREO-B, "/V" for VEX. Column R indicates the range of 
solar radii where the measurements were obtained or techniques were applied. For the 
fitting methods, this is the range of the distances Rfp and Rhm corresponding to the first 
and last e(t) measurement of the ICME track in the Jmap. Angle I is the longitude of the 
source region (flare peak), the ICME direction (FPF/HMF) and the spacecraft positions 
(STEREO-B I VEX) with respect to the STEREO- A spacecraft (I > for solar west). V 
is the CME/ICME speed, and t / t a are the launch and arrival times, respectively. 



technique 


Instruments 


R Q 


I, deg 


V, km s 1 


to, UT 


ta, UT 


flare peak 


GOES 


1 


-89.5 




Feb 13 05:46 




tracking 


COR1A 


2.5 


-91 a 


350 a 


Feb 13 06:20 


Feb 18 03:57 b 


FPF/B 


HI1/2A 


14.1-203.7 


-79 


280 


Feb 13 01:30 


19 Feb 06:12 


HMF/B 


HI1/2A 


14.7-225.7 


-107 


322 / 310 c 


Feb 13 03:20 


18 Feb 18:00 


STB 


IMPACT/PLASTIC 


215.6 


-91 


362 / 348 d 




18 Feb 10:00 


FPF/V 


HI1/2A 


14.1-203.7 


-79 


280 


Feb 13 01:30 


17 Feb 12:01 


HMF/V 


HI1/2A 


14.7-225.7 


-107 


322 


Feb 13 03:20 


18 Feb 02:36 


VEX 


MAG 


154.5 


-68.9 






17 Feb 22:00 



a The COR1A speed was measured in the plane-of-sky by Miklenic et al. (2011, submitted to Solar Physics). 
fc Ballistic projection of the CMEs arrival at STEREO-B with the COR1 measurements. 
c The first speed is for the apex of the HM circle (Vrm) and the second for the flank part hitting STEREO-B 
(Vhmi)- 

d The in situ V are means over the ICME sheath / flux rope regions. 
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Fig. 1. — Illustration of the different elongation-to-distance conversion methods. From the 
top left, an observer looks towards a circular ICME front. The green (PP), red (FP) and blue 
(HM) dots show where the observer would calculate the radial distance of the ICME from 
the Sun using the indicated methods. In this case, the observer used the same direction <f> for 
the FP and HM conversion equations, which leads to Rfp > Rhm for the same elongation 
value. 
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Fig. 2.— Evolution of the ICME in STEREO-A HI1 (top 2 images) and HI2 (bottom 2 
images). The ICME is delimited in elongation e and elevation 5 by dashed and dotted lines, 
respectively. These coordinates represent longitude (e) and latitude (5) in a Helioprojective 
Cartesian (HPC) coordinate system. In the left two images the ecliptic (in the plane-of-sky 
of STEREO-A) is plotted as a dashed- dot- dot line. The positions of STEREO-B (triangle) 
and Venus (cross) are indicated in the bottom two images. In the right two images the 
upper, central, and lower position angles of the ICME are indicated, again for the plane- 
of-sky of STEREO-A, giving a directly measured angular width of AWo = 23° (HI1) and 
AWn = 19° (HI2). The angular widths AWl in all images are based on the formula by 



Lynch et al. 



( 120101 ). and are quite similar to the AWd values. The small angular width and 



the morphology are consistent with a CME viewed edge-on, so its axis of symmetry points 

nut nf thp Tinner nlanp This fitnirp is alsn availablp a.s an mnpp - animation in thp plpptrnnir 
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Fig. 3. — (a) Proton density observed by STEREO-B/ 'PLASTIC. Enhanced densities (up 
to 25 protons cm -3 ) in front of the flux rope are delimited by two vertical lines, with the 
left one setting the arrival time, (b) Jmap along the ecliptic plane with the elongation of 
STEREO-B marked as a white horizontal line. It is seen that a track, starting on Feb 13, 
reaches the elongation of STEREO-B and Venus on early Feb 18, which matches very well 
the single density jump observed by STEREO-B. (c) Best HMF / FPF fits for the ICME 
leading edge. Both functions are able to reproduce the observed function e{t). (d) Residuals 
between the fitting functions and the observations. The time development of e(tj) — £fu{ti) 
is similar for both techniques and the summarized residue a 2 for HMF is half the value for 
FPF. 



-30 - 



STEREO-B IMPACT / PLASTIC 17-20 February 2009 A t=l min 




10 

I 

PQ 

-10 
500 ■ 

| 400 

> 300 : 



— R 




— T 




N 





Fixed Phi 



arrival 



Harmonic Mean 



<- flux rppe -> I 




o.oi 

10 

1 

0.1 
0.01 



^ — ^-f~^l 1 1 
1 1 1 1 1 1 




1 


i i i i i 
i i i i i 














i 






00 06 12 18 00 06 12 18 00 06 12 18 00 06 12 IS 



02/17 



02/18 



02/19 



02/20 



Fig. 4. — Magnetic field and plasma data for STEREO-B. The first solid line marks the 
approximate arrival time of the density front. The 2nd and 3rd solid lines indicate the flux 
rope interval. The arrival times from the FPF and HMF methods are given by red and blue 
vertical dashed lines, respectively, and the fit velocities Vfp and Vhmi are given by red and 
blue horizontal lines, which indicate the extent of the arrival time error ±54 min (FPF) 
and ±20 min, which are derived from the complete ICME track. From top to bottom: 
magnetic field magnitude and components in RTN coordinates (R pointing radially away 
from the Sun, T is the cross product of the solar rotation axis and R, N completes the right- 
handed triad), proton bulk velocity, proton number density, proton temperature (black) and 
expected temperature (red), proton (3 and the total, magnetic and plasma pressure. 
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VENUS EXPRESS Magnetometer 15-20 February 2009 
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Fig. 5. — Venus Express magnetometer observations on 2009 February 15-20. Originally, the 
magnetic field components are given in Venus solar orbital coordinates, where X is directed 
towards the Sun, Y is opposite to the planets orbital velocity (?» solar east), and Z points 
towards north completing the right handed triad. To make a comparison easy with Figure HJ 
we plotted here —B x and — B y components so that they approximately match the orientation 
of the RTN coordinate system (though for the position of Venus) used for the STEREO-B 
magnetic field components. The arrival times at Venus are marked again by red (FPF) and 
blue (HMF) dashed lines. An ICME interval of higher than average magnetic field strengths 
and smoothly rotating field from Feb 17 22:00 UT - Feb 18 14:00 UT is highlighted (shaded). 
For comparison, the arrival time at STEREO-B is indicated by a solid vertical line. 
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ICME 13-18 Feb 2009 
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Fig. 6. — Illustration of the ICME directions inferred for the 2009 February 13-18 ICME 
event using the FPF and HMF methods, looking from north down onto the ecliptic plane. 
The positions of the STEREO spacecraft, Sun and Earth are indicated. The black (green) 
circles gives the shape of the Harmonic Mean front at the arrival time of STEREO-B (Venus). 
Because the HMF method assumes a circular front, the calculated arrival times at both 
STEREO-B and Venus matches better the observed arrival times. Two black arrows show 
the rotation of the magnetic field inside the flux rope at STEREO-B, its radial size indicated 
by a green cylinder. 
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ICME arrival time as function of track length 
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Fig. 7. — Accuracy of possible predictions for the ICME arrival time (upper panel) and 
direction (lower panel) with respect to the arrival time and location of STEREO-B, using 
different lengths for the tracks extracted from the Jmap, given by their maximum elongation 
value max(e). Top: The difference of the arrival times t a calculated for the FPF (red) and 
HMF (blue) methods to the observed arrival time £g at STEREO-B. Indicated are also the 
outer boundary for the HI1 FoV (dashed vertical line) and the inner boundary of the HI2 
FoV (solid vertical line). Bottom: Differences between the calculated propagation directions 
and the longitude of STEREO-B (As = — Ib), as function of max(e). 
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A. Alternative version of the inverted Harmonic Mean formula 

The aim of this section is to explicitly show how to solve Equation ([3]) for e(t). This is 
the conversion formula from elongation to distance if one assumes that the observer looks 
along a tangent to a circle which is always attached to the Sun at one end ("Harmonic 
Mean" assumption, see Figures 1 and[8k). The radial di stance of the circle 's apex Rhm, 



the point with largest distance from the Sun, is given by 



Lugaz et al. 



(I2009h : 



2ri sine(t) 

riHM{t - — — . , (At J 

1 + sm(e(t) + <p) 

with d the radial distance of the observer from the Sun, e(t) the observed elongation as 
a function of time and the constant ICME propagation direction with respect to the 
observer (0 > means solar west). After some basic algebra, Equation (IA10 can be written 

as 

( ?f° . - cos (pHM J sin(e(t)) - sin(0# M ) cos(e(t)) = 1, (A2) 

V VHMt J 

where we have already replaced distance with speed: Rhm = Vnut- For better readability, 
we define 

2d a 

a = — — — - cos((Phm); b = sm((p HM ), (A3) 



VHMt 



which results in 



asin(e(t)) - 6cos(e(t)) = 1. (A4) 
This equation can be brought into quadratic form by re-arranging it to 

ay/1 -cos 2 (e(t)) = 1 + 6 cos(e(t)) (A5) 

and squaring both sides. Solving the resulting quadratic equation in cos (e(t)) leads to our 
form of the Harmonic Mean fitting equation as 

e(t) = arccos( b + a ^ a + b 1), (A6) 
a 2 + b z 
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We conclude this section with some notes useful when developing a code that fits the 
inverse equations. While the arctan function in the FPF Equation (J2J) is symmetric with 
respect to positive and negative values of the direction 0, the HMF equation (IA6I) will 
always produce positive values for e(t) because of the arccos function, regardless of the 
sign of the direction <Pfp- Elongations provided by STEREO-A/HI, which looks to solar 
east, are negative by definition, so these e < directly translate to (ppp < (or solar east) 
when applying the FPF method. For STEREO-B , automatically e > and thus 4>fp > 0, 
consistent with our definition that 4>fp > corresponds to solar west. In contrast, when 
applying the HMF method to STEREO-A, always positive elongations have to be fitted 
and the sign of the direction 4>hm has to be additionally corrected. The reason for this 
asymmetry is that the observer is also able to measure an ICME front with HMF even if 
the apex of the HM circle is on the other side of where the Helio spheric imag er is looking, 



i.e. when the apex is west of STEREO-A or east of STEREO-B (ILugaz 



2010). 



B. Arrival time calculation for circular ICME fronts 



For calculating the arrival time of a circular ICME front at a given 



the heliosphere, as assumed by t 



tangent-to-a-sphere (ILugaz et al 



le Ha rmonic Mean (HM) fitting (ILugaa 120101 ) or the 



posit ion in 



20101 ) methods, it is necessary to derive a formula which 



takes into account that different parts of the front will move with different speeds. As a 
consequence, the ICME arrival time will be a function not only of the radial distance of the 
in situ spacecraft / (or planet) from the Sun but also of its separation angle A to the HM 
apex, the point with the largest distance from the Sun. 

If the ICME in question is not heading exactly towards the in situ observing spacecraft 
/, i.e. the ICME direction <Phm as t ne result of the Harmonic Mean fitting method is not 
the same as the longitudinal separation between the observer O and /, there will be a time 
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Sun 



(b) 

Arrival time correction for circular CME front 




A [degrees] 

Fig. 8. — (a) A sketch showing the definitions and geometry for deriving the arrival time 
correction for the HMF and TAS methods (see text). The solid circle is the circular front 
assumed for the ICME leading edge and the dashed line depicts 1 AU. (b) The arrival time 
correction - Equation ( 1B5I) - is plotted as a function of the angular separation between the 
ICME apex and the in situ observing spacecraft /, which is situated at a radial distance 
from the Sun of di = 1 AU. 
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delay t c > when the ICME hits I compared to when the ICME apex is at the same radial 
distance of / from the Sun. Figure [Sh^ visualizes this geometry. We will first derive this 
formula as a "correction" t c compared to the ICME apex to assess its significance, and will 
then derive the formula useful for direct calculation of the arrival time t a HM- The arrival 
time correction t c is a function of the angular separation A = 4>hm — h between 4>hm and 
the heliospheric longitude of I (called 1+), the velocity Vhm obtained by the fitting procedure 
and the radial distance di of the in situ spacecraft from the Sun: t c (A, Vhm, di). We use 
all longitudinal angles as such that westward is defined as positive, and the condition 
A G [— 7r/2... + 7r/2] is necessary for the spacecraft to be hit by the expanding circular front. 
Additionally, c is defined as the distance the HM apex travels before the HM circle hits the 
spacecraft I. Applying the law of sines to the triangle indicated in Figure [SK, we find 

(di + c)/2 di 



sin(A) sin(7r - 2A) ' 

Simple regrouping leads to 



(Bl) 



c = di( — rr - 1 . (B2) 
Vcos(A) J 

To convert this distance to a time, we assume that the ICME continues to travel with the 
constant velocity obtained from the fitting procedure (Vhm) and define t c > as the arrival 
time correction, 

U = c/V„ M (B3) 

t c = ' UC0S ' 1{A> - L > (B4) 
Vhm 

This can be rewritten as the simple formula 

t c h = 41555 x ,L U 1 _ >- B5 

V HM [km s l \ 

with the result t c in hours and Vhm an d di given in convenient units. Figure |Hb shows the 
behavior of t c (A) for di = 1 AU and various Vhm- Obviously, for faster CMEs the arrival 
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time correction becomes less, and for a Vhm = 800 km s -1 CME, the arrival time correction 
is less than 10 hours for A < 35 degrees. It might be surprising that Equation (IB5I) depends 
on di - this is simply because the circle assumed by HM is always attached to the Sun at 
one end and thus expands while the apex propagates outward. 

In summary, to calculate the arrival time of an ICME at a given heliospheric distance 
di and at a longitudinal separation A with the HMF method, a correction t c must be added 
to the arrival time of the ICME apex: 



taHM — toHM + TJ~ ^ 0^) 

VHM 

with toHM being the launch time (defined here by eitoHu) = 0). This reduces with 
Equation flB4j) to simply 



taHM — toHM + 77 77TT- (B7) 

Vhm cos(A) 

This is the final equation for calculating the arrival time of an ICME with the HMF 
method. From the last equation, the elementary result follows that the speed of the ICME 
flank at a given angle A to the ICME apex is just reduced by a factor cos(A) for a circular 
geometry. When comparing the speed Vhm as an outcome of the HMF method to in situ 
measurements of the solar wind bulk speed, the speed Vhmi = Vhm cos(A) has to be used. 
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